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SUMMARY 


Determining the three-^dimensional magnetostatic field in current-induced 
situations has usually involved vector potentials, which can lead to excessive 
computational times. A recent paper shows how such magnetic fields may be 
determined using scalar potentials. The present paper shows how the heat 
transfer capability of NASTRAN. Level 17 has been modified to take advantage of 
the new method. 

INTRODUCTION 


All classical electromagnetic phenomena are governed by 
equations: 

the four Maxwell 

VXH- J + ll 

(1) 

V xE + 1“ = 0 
d tZ 

(2) 

V • D = p 

(3) 

V • B = 0 

(4) 


where the vector quantities are defined as follows; 

H Magnetic field strength or intensity 
B Magnetic Induction or flux density 
J Current density 
E Electric field strength 

D Electric displacement 

and the scalar quantities are defined as follows; 

p Charge density 

t Time 


There is also a constitutive relation between B and H, given by 


(5) 


B = pH 
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where p is the magnetic permeability. 

Electromagnetic problems are often solved by Introducing and solving for 
the magnetic vector potential A, where 

B = VxA (6) 

Spreeuw and Reefman (tef. 1) used this nsethod with NASTRAN in solving for 
harmonically oscillating electromagnetic fields in the presence of conductors 
carrying alternating currents. However, in order to use the existing 
structural and heat transfer capabilities in NASTRAN, simplifying assumptions 
had to be made. In particular, the magnetic vector potential A and the source 
current densities J were allowed to have components in only one direction, and 
those components were invariant in that direction. These assumptions 
effectively reduce the problem to one of solving for a scalar potential, which 
can be handled by NASTRAN' s heat transfer analyses. 

In the same paper, Spreeuw and Reefman also considered a problem in which 
A was not unidirectional and were able to use NASTRAN 's structural analysis 
capability onlv because the governing equations uncoupled for the components of 

A. 


Another problem with this formulation is the requirement that 

V • A = -£p|| 

where e is the electric permittivity. 

Spreeuw and Reefman used a separate post-processor to handle this condition. 

Frye and Kasper (ref. 2), in solving magnetostatic problems using the vector 
potential, use a Lagrange multiplier method (similar to multi-point constraints) 
to satisfy constraint (7) . They also point out that special boundary conditions 
are required at boundaries where the permeability p changes. 

In reference 3, Zienkiewicz, Lyness, and Owen have developed a method for 
solving general, three-dimensional magnetostatic problems using scalar 
potentials, so that standard heat transfer analyses may be used and constraint 
equation (7) is not required. They also indicate that special boundary condi- 
tions, such as those mentioned in reference 2, are not needed. 

The present paper shows how this new method has been implemented in 
NASTRAN Level 17. 


BASIC EQUATIONS AND ASSUMPTIONS 


The problem to be solved is the determination of the magnetostatic field 
due to a body placed in an existing magnetic field produced, for example, by 
direct current-carrving loops. The materials are assumed to be linear, but 
may be anisotropic. The governing equations are; 

7 xH = J (8) 
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B «= 

yH 


(9) 


y • B 

n 

0 

(10) 

Zienkiewicz separates U into two 

parts , 





H « H 


H 

(11) 


c 


m 

Hj, is the field In a homogeneous 

region 

due to current J, satisfies 



y xH 

c 


J 

(12) 


and is computed using the Biot-Savart lav/. is the unknovm magnetic field 
strength and satisfies 


VxH =0 (13) 

m 

so that Hjjj = (14) 

and V • y V<|. + V . p = 0 (14a) 

where ^ is the scalar potential. Zienkiewicz then uses standard variational 
principles, with equations (9) and (10), to arrive at the standard finitt 
element form 


K(j) = F 

where K is the "stiffness" matrix, 

F is the "load" vector, and 



f = J (VN II dV 
1 „ i c 


(15) 


(16) 

(17) 


where 

til 

N. is the finite element shape function for the i grid point, and 
i 

V is the volume of the finite element. 

The formulation (16) for k^j is exactly that of the standard heat transfer 
conductivity matrix with magnetic permeability y playing the role of thermal 
conductivity. The formulation (17) for f^, however, is not a standard heat 
transfer quantity and must be computed either in a separate program and input 
to NASTRAN or in a new NASTRAN capability. Also, note that fi is element- 
dependent, as evidenced by the shape function gradient in equation (17) . 
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Equation (15) is solved for <{i subject to standard natural or forced 
boundary conditions, Hjn can then be computed from equation (14), and the final 
results can be obtained using aquations (11) and (9) . 


NASTRAN 1^^PLEMENTATI0N 


To solve magnetostatic problems with NASTRAN Level 17 using the methods of 
the previous section, we select rigid format 1, HEAT approach. However, we 
have modified the program to 

1) compute the li^ field due to circular, direct current-carrying loops; 

2) accept a specified Hq field; 

3) compute fj^ (equation (17)) for the axisymraetric solid ring elements 
TRAPRG and TRIARG; 

4) perform the addition specified in equation (11) , where H[n is a 
transformation of the "temperature" gradients; and 

5) output B (equation (9)) for subsequent NASTRAN plotting. 

The implementation thus far has been limited to solid axisymmetric 
problems using TRAPRG and TRIARG finite elements and is running on the DTNSRDC 
GDG 6000 computers . 


NEW BULK DATA CARDS 


Two new bulk data cards have been introduced into the program for 
computing fields. They are CEMLOOP, for computing the Hj, fields due to 
circular current loops, and SPCFLD, for specifying at selected grid points. 
(See figures 1 and 2 for detailed descriptions of these cards.) 


MODIFIED NASTRAN ROUTINES 


Nineteen existing NASTRAN routines have been modified to accommodate the 
new capability. The routines and the nature of the modifications are as 
follows; 

Routine Reason for Change 

IFF , IFS4P , IFX1BD-IFX7BD Recognize and check new bulk data cards 

CEMLOOP and SPCFLD. 

LD21 Restart with CEMLOOP and SPCFLD in Static Heat 

Transfer Analysis. 

GP3A,GP3BD Recognize CEMLOOP and SPCFLD as "heat transfer" 

load specifications and place on HSLT (Heat 
Static Load Table). 
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Routine 

SSGSLT 


EXTERN 


XBSUD, X^ffLlJD 

XSE^aA 

OFPIA 


OFP2ZZZ 


Reason l:oi.‘ ChaiiRe 

Retrieve CEMLOOP end SPCPLD specifications from 
IISLT. CEMLOOP cards are copied directly to 
data block NEWSLT for later processing. All 
SPCFLD specifications are combined into one 
vector giving total Hj, at all grid points. 

This vector is placed on NEWSLT. 

Call new routine EANDM, which retrieves CEMLOOP 
and SPCFLD specifications from NEWSLT and 
computes f as given in equation (17) , 

Specify table updates so that NASTRAN will 
recognize new functional module EMFLD. 

Call new functional module EMFLD. 

Specify new headings to output new data block 
HOEHl giving magnetic field strength and 
induction. 

Call for new headings from OFPIA when HOElll is 
recognized , 


NEW ROUTINES 


Three new routines have been developed. Subroutine EANDM reads into open 
core CEMLOOP and SPCFLD data from data block NEWSLT, reads information for an 
element from the Element Summary Table (EST) , and calls an element-dependent 
routine to compute the load os given in equation (17) . 

Subroutine EMRING is called by EANDM aiid computes the load due to CEMLOOP 
and SPCFLD specifications for solid axisymraetric trapezoidal (TRAPRG) and 
triangular (TRIARG) rings. We assume that is constant over an element. 
Therefore, for TRIARG elements, is computed at the centroid due to all 
CEMLOOP 's using an elliptic integral formulation. If SPCFLD 's are given, llj, is 
computed to be the average of the Uq's at the three vertices. Each TRAPRG 
element is divided into four overlapping triangles, and each triangle is treated 
Of.’ a TRIARG element. Once Hj. has been computed for a triangle, equation (17) 
is used to compute the load at each grid point forming the triangle. Subroutine 
EMRING also outputs to Fortraia file II certain elemev\t information, including 
H(-, for later processing by functional module EMFLD. This is a temporary 
method for passing information from EMRING to EbfFLD. The normal method, of 
course, is to use data blocks. However, subroutine EMRING resides in functional 
module SSGl, and adding a nev^ output data block to that existing module would 
require a change to every rigid format in the program. These changes will be 
made at a later time. 

Subroutine EMFLD, which is also a new functional module, computes the 
magnetic field strength and induction, according to equations (11) and (9), for 
each finite element in the model. EMFLD retrieves from data block HOEFl the 
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"temperature" gradient for an element. Since the gradient U,„ waa computed in 
an element coordinate system, EMFLD transforms it into basic coordinates. Then 
Fortran file 11 is searched to match the element identification number, and, on 
a match, li^ for the element is retrieved, added to the transformed Hi,|, and put 
out to data block HOEHl. Also computed and output to HOEHl is the magnetic 
induction B. HOEHl is later output using normal Output Pile Processor (OFF) 
execution, EMFLD also computes and outputs other information for plotting 
purposes, as explained in the next section. 


PLOTTING MAGNETOSTATIC RESULTS 


Normal NASTRAN plot processing allows for deformed plots based on grid 
point displacements or contour plots of stresses. In the present analysis, 
however, the "displacements" (the scalar potentials) are of little use by 
themselves. The "stresses" (H or B fields coming from EMFLD) are more useful, 
but what we would really like to see for "nice" plots are the lines of 
induction. The lines of induction indicate the direction of the magnetic 
induction B, and the number of lines per unit area indicates the magnitude of 
B. Wliile we do not presently plot these lines of induction, we do plot the 
actual induction, magnitude and direction, for each element. Therefore, 
functional module EMFLD outputs other quantities as follows. For each element, 
two coincident grid points are created at the centroid of the element, and the 
corresponding GRID cards are punched. Also punched is a PLOTEL card connecting 
the two grid points. (The length of the PLOTEL element is zero.) Then a 
"displacement" vector is created by assigning zero values to each of the 
original grid points in the model and assigning the magnetic induction value 
for an element to each of the two coincident grid points created for the 
element. (The "displacement" vector uses six degrees of freedom per grid point 
since B is a vector, not a scalar.) This vector is packed in EMFLD and output 
in DMAP using module OUTPUTl. On a subsequent NASTRAN run, the new GRID and 
PLOTEL cards are merged with the original data, the "displacement" vector is 
retrieved using DMAP module INPUTTl as data block UGV, and a deformed plot is 
requested with the VECTOR R option. This NASTRAN run is performed with rigid 
format 1, DISP approach, and ALTERs are used to execute only those modules 
required for deformed plots. The plots show the original structure as an 
underlay, and a vector is drawn at each element centroid indicating the magni- 
tude and direction of B in that element. 


SAMPLE PROBLEMS 


At the time that this paper was being prepared, the only axisymmetric 
problems run with NASTRAN for which analytical results were readily available 
were problems with uniform permeability. The comparison between the NASTRAN 
and analytical results was very good. Although problems with nonuniform 
permeability have been run with NASTRAN and "reasonable-looking" results have 
been obtained, analytical results, required for verification, are still 
forthcoming. 
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Input Data Card CEMLOQP Circular Current Loop 


Description; Defines a circular current loop in magnetic field problems- 
Format and Example; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CEMLOOP 

SID 

J 

AXI 

XI 

Y1 

Z1 

X2 

Y2 

+A 

+A 

Z2 

xc 

YC 

zc 

CID 





CEMLOOP 

3 

2.5 

1 

1 

5.2 

0. 

2.25 





Field 

SID 

J 

AXI 

X1,Y1,Z1 

X2,Y2,Z2 

XC,YC,ZC 

CID 


Contents 

Load set identification number (integer > 0) , 

Current through loop (units of positive charge/sec) (real i 0) . 

= 0 nonaxisymmetric problem (not yet implemented) 

= 1 axisymmetric problem; TRAFRG and TRIARG elements are 
implied (integer) , 

Coordinates of two points through which the loop passes 
(given in coordinate system CID) (real) . 

Coordinates of center of loop (given in coordinate system 
CID) (real) . 

Coordinate system identification number (Integer t 0) . 


Remarks : 


1. Load sets must be selected in the Case Control Deck (LOAD=yiD) to 
be used by NASTRAN . 

2. If AXI=1, Y1 must be 0. or blank, and all data fields after Z1 must 
be 0. or blank. (Continuation card need not be present.) 


Figure 1. Bulk Data Description of CEMLOOP 
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Input Data Card SPCFLD Specified Magnetic Field 


Description; Specifies magnetic field at selected grid points. 
Format and Example; 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

SPCFLD 

SID 

HCX 

HCY 

HCZ 

G1 

G2 

G3 

G4 


SPCFLD 

18 

12.25 

0. 

62. 

8 

1/ 

103 

1 



or 


SPCFLD 

SID 

HCX 

HCY 

HCZ 

GlDl 

THRU 

GID2 



SPCFLD 

18 

12.25 

0. 

62. 

9 

THRU 

27 




or 


SPCFLD 

SID 

HCX 

HCY 

HCZ 

“1 





SPCFLD 

18 

12.25 

0. 

62. 

“1 






Field C ontents 

SID Load set identification number (integer > 0) . 

HCX,11CY,HCZ Components of specified field (real). 

GijGIDi Grid point identification numbers (integer > 0). 

-L Implies the field applies at all grid points. 

Remarks ; 

1. Load sets must be selected in the Case Control Deck (L0AD=SID) to 
be used by NASTRAN . 

2. All grid points referenced by ClDl THRU GID2 must exist. 


Figure 2. Bulk Data Description of SPCFLD 
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